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Abstract 



The efficiency of different types of Brownian motors is calculated analytically and nu- 
merically. We find that motors based on flashing ratchets present a low efficiency and an 
unavoidable entropy production. On the other hand, a certain class of motors based on adi- 
abatically changing potentials, named reversible ratchets, exhibit a higher efficiency and the 
entropy production can be arbitrarily reduced. 

In the last years there has been an increasing interest in the so-called "ratchets" or Brownian 
motors ||, ||, ||, Q |j| These systems consist of Brownian particles moving in asymmetric 
potentials, such as the one depicted in fig. [l] (left), and subject to a source of non-equilibrium, 
like external fluctuations or temperature gradients. As a consequence of these two ingredients 
— asymmetric potentials and non-equilibrium — , a flow of particles can be induced. 

Most of the cited papers consider systems where the Brownian particles do not gain energy 
in a systematic way. Although these systems are called "Brownian or molecular motors" , they 
do not convert heat into work, nor induce any energy conversion. Feynman in his Lectures 
already understood that, in order to have an engine out of a ratchet, it is necessary to use its 
systematic motion to store potential energy. This can be achieved if the ratchet lifts a load. Then 
the ratchet becomes a thermal engine and Feynman estimated its efficiency (although following 
assumptions which have been revealed to contain some inconsistencies ||). Recently, Sekimoto 
H has generalized this procedure, defining efficiency for a wide class of ratchets. Jiilicher et 



al H have also discussed the efficiency of molecular motors and Sokolov and Blumen 1 1 have 
calculated the efficiency of a deterministically flashing ratchet in contact with thermal baths at 
different temperatures. A general conclusion is that these motors are intrinsically irreversible, even 
in the quasistatic limit || ^, [l^] . 

On the other hand, it has been recently introduced |llf| a class of deterministically driven 
ratchets where the entropy production vanishes in the quasistatic limit, i.e., reversible ratchets. 
The aim of this letter is to explore the differences, regarding efficiency, between randomly flashing 
ratchets and both reversible and irreversible deterministically driven ratchets. 



Randomly flashing ratchets 

Consider two species of Brownian particles, say A and B, moving in the interval [0, L] with periodic 
boundary conditions. Particles of type A feel a potential Va(x), whereas particles B feel Vb{x). 
Besides, there is a continuous transfer of particles, A ^ B, which accounts for non-equilibrium 
fluctuations. This picture is equivalent to that of a single Brownian particle in a randomly switching 
potential ||. 
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Figure 1: Left: Asymmetric sawtooth potential of the ratchets presented in refs. [[ij |[ || [|, || . In 
this letter we consider two types of ratchets: a) one where the potential is randomly switched on 
and off; and b) one where the potential is deterministically modulated. 

Right: Efficiency and current of the ratchet where the the potential on the left is randomly switched 
on and off (case a), as a function of the maximum height V of the potential. The reaction rates 
are uja = 1-08 and uj b = 81.8, a = 1/11, and the external force is F = 4.145. 



In ref. ||, it was proved that a flow towards a given direction, say, to the right, occurs for some 
asymmetric potentials Va and Vb ■ If we add a load or force F opposite to the flow, the evolution 
equation for the probability density of particles A, pa{x), and particles B, p B (x), is: 

d t pA{x,t) = -d x J A pA{x,t) - oj\p A (x,t) - p B (x,t)] 

d t pB(x,t) = -d x J B pB(x,t) +u>[p A (x,t) - p B (x,t)] (1) 

where Ji — —Vl(x)—F — d x is the current operator, the prime indicates derivative with respect to 
x, and d is the rate of the reaction A ^ B. We have taken units of energy, length and time such 
that the temperature is k B T — 1, the length of the interval is L = 1, and the diffusion coefficient 
is D = 1. 

The flow of particles in the stationary regime is J = JaPa\x) + JbP s b {x), where p A l B (x) are 
the stationary solutions of eq. ([!]) . This flow J is a decreasing function of the external force F and 
becomes negative if F is stronger than a balancing force, Fj, a i. Therefore, if < F < F^ai, particles 
move against the force and, consequently, gain potential energy in a systematic way. The potential 
energy gain or output energy per unit of time is E out — JF, which vanishes both for F = and 

F = Fbal- 

On the other side, switching on and off the potential requires some energy. In our two-species 
picture, the reaction A ^ B does not conserve energy since V A (x) ^ V B (x). Therefore, in each 
reaction A — > B, occurring at a point x, an energy V B (x) — V A (x) is transferred to the system 
(or withdrawn, if the sign is negative). Similarly, an energy Va(x) — V B (x) is transferred to the 
system in each reaction B — > A occurring at x. In the stationary regime, the average number of 
such reactions per unit of time is, respectively, ujp A t (x) and u>p B (x). Therefore, the input energy 
per unit of time is |(| |j| : 

E in =u f dx [V B (x) - V A (x)] [(4(x) - p s i{x)] . (2) 
Jo 

Finally, the efficiency can be defined as || ^| : 

V = -p— • (3) 

This efficiency can be calculated analytically for the system given by eq. (||) with piecewise 
potentials. We have performed an exhaustive study for the particular setting V B {x) =0 and Va(x) 
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equal to the potential depicted in fig. |T] (left): 



„ , , f Vx/a if x < a . 

V ^ = { V(l-x)/(l-a) Xx>a ^ 

with a = 1/11. For this system, the maximum efficiency is rj max = 3.29%, which is reached for 
V = 22, F = 3, and lu = 63. This efficiency can be improved using different reaction rates: u>a for 
A -> B and w B for B -> A. In this case, i] max = 5.315% with V = 16.7, F = 4.145, uj a = 1.08, 
and = 81.8. Observe that, with these values for uja,b, the particle stays much longer within 
the potential Va(x) than within Vb{x). 

We have plotted in fig. |l| (right) the efficiency and the flow of particles as a function of V 
setting the rest of parameters equal to these optimal values. Two are the messages from this 
figure. Firstly, the maximization of the efficiency is a new criterion to define optimal Brownian 
motors and this criterion is, in some cases, less trivial than that of maximizing the flow. Secondly 
the randomly flashing ratchet under study has a rather low efficiency. As we have mentioned 
before, the heat dissipation per unit of time is Ei n — E out . Consequently the increase of entropy 
of the thermal bath, per unit of time, is Ei n — E out , since fceT = 1. On the other hand, in the 
stationary regime there is no change of entropy in the system nor in the external agent which 
provides the non-equilibrium fluctuations^. Therefore, the net entropy production per unit of time 
is Ei n — E out . If this entropy production would vanish, i.e., if the system could work in a reversible 
way, it should reach a 100% efficiency. On the contrary, the efficiency is below 10% and we can 
conclude that the motor based on the randomly flashing ratchet is very inefficient. 

One could think that the efficiency would increase in limiting situations where the system is 
close to equilibrium, such us u> — > and/or V A ~ Vb — * 0. However, a perturbative analysis of 
eq. ([j]) shows that r\ — > in both limits. In the first case, to — > 0, from eq. (|l|) one can easily 
find that J is of order oj, so is Fbai- Therefore, E out , in the interval < F < Fb a i, is of order u 2 , 
whereas one can prove that Ei n is of order to, giving a zero efficiency in this limit. In the second 
case, AV(x) = V A (x) — Vb(x) — ► 0, the input energy E in is of order AV 2 . However, surprisingly 
enough, J is of order AV 2 and so is F^ai, yielding E out of order AV 4 and, again, a vanishing- 
efficiency We conclude that the flashing motor is intrinsically irreversible, as it has been pointed 
out for related models in refs. §, §, [| flOf. 



Deterministically driven ratchets 

A different strategy to reduce the production of entropy consists of considering Brownian particles 
in a potential which changes deterministically in time. If the potential is changed very slowly, the 
system evolves close to equilibrium and the entropy production is low. From now on, we focus 
our attention on Brownian particles in a spatially periodic potential V(x; R(i)) depending on a 
set of parameters collected in a vector R which changes in time p| . The parameters are changed 
periodically in time with period T, i.e., R(0) = R(T). 

As in ref . [pi , we have to modify our definition of efficiency. Firstly, we deal with energy transfer 
per cycle [0, T] instead per unit of time. Secondly, the input energy or work done to the system in 
a cycle, as a consequence of the change of the parameters R(i), is: 

f T , dV(x;R(t)) . , 
E m = dt g t POM)- (5) 

The probability density p(x, t) verifies the Smoluchowski equation: 

dtpjx, t) = -d x J R{t) p(x, t) (6) 

1 A physical realization of this external agent is a third species of particles, say C, feeling a potential Vc(x) = 
Vb( x ) ~ Va( x ) an d participating in the reaction as A + C ^ B. If the temperature of C particles is the same as B 
and A particles, then detailed balance holds and there is no flow of particles. However, if the temperature of the C 
particles is infinity, we recover the flashing ratchet discussed in the text. Therefore, this randomly flashing ratchet 
can be considered as a therm al e ngine in contact with two thermal baths, one at T = 1/kg and the other one at 
infinite temperature (see also flOl for an interpretation of the deterministically flashing ratchet as a thermal engine) . 
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where J7r = — V'(x] R) — F — d x is the current operator corresponding to the potential V(x; R). 
As before, the output energy is the current times the force F, but now the current is not stationary 
and we have to integrate along the process: 

E out = [ dt FJ m) p{x,t) = F(f>. (7) 
Jo 

where <fi is the integrated flow. 

With the above expressions, the efficiency of the system, r\ = E out /Ei n , can be found analyti- 
cally for T large and weak external forces, where it is expected to be high. For the integrated flow 
one finds <p = <po — pFT, where p is the average mobility of the system: 

_ _ 1_ f T dt 

"~tJ Z + (R(i))Z_(R(i)) ( > 

and 4>q is the integrated flow for F — [[llj : 

00= &dH- [ dx I dx'p + (x;R.)V R p.(x';R), (9) 
J Jo Jo 

with 

p ±V(x;R) rl 

P±(x;R) = ~z^jf£j~'> Z ±( R ) = J dx e ±v ^ R \ 

In eq. (^|) the contour integral runs over the closed path {R(t) : t G [0, T}} in the space of 
parameters of the potential. The term proportional to T in the integrated flow, cj> — 4>o — pFT, 
arises because the force F induces a non-zero current which is present along the whole process. 
As a consequence, the balancing force is F^ a i = (j} /pT, and, in order to design a high efficiency 
motor, it is necessary to take simultaneously the adiabatic limit T — > 00 and the limit F — ► with 
FT finite. Notice also that the above expressions are useless if 0o = 0. In a previous paper pj] ], we 
have discussed the conditions for </>o to be different from zero and called reversible ratchets those 
systems where <fio ^ 0. From now on, we restrict our analytical calculations to reversible ratchets, 
although we present below numerical results for an irreversible ratchet. 

The input energy for weak force F and large T is Ei n = 4>qF + b/T, with 

b = - f dt Z_(R(t))Z+(R(t)) ( f dx f dx' p + (x;R(t))[dtp^(x';R(t))} 
Jo { Uo Jo 

+ [ dx f dx' [ dx" [d t p-{x;TL{t))] p + (x';R(t))[d t p-(x";R(t))]\ (10) 
Jo Jo Jo J 

which is a positive quantity. Combining the above expressions, one finds for the efficiency: 

_ F((/)o - pFT) _ cf)Q a - pa 2 
V ~ <poF + b/T 4>oa + b [ ' 

where a = FT. This expression is exact in the limit T — > 00, F — > 0. Notice that, even for large 
T, the irreversible contribution, b/T, to Ei n is of the same order as (j>oF. 

In a given system, i.e., for a set of parameters <fio, p and b, the maximum efficiency is reached 
for a — (b/cj)o)[y/l + (jtfjjflb) — 1] and its value is given by 



Vmax = 1-2 y/z(l + Z) - Z (12) 



with z = bp/4>Q. Eq. (|l2j) clearly shows how the term b in the denominator of eq. ( |ll| ) prevents 
the system to reach an efficiency equal to one. Fortunately, as we will see below in a concrete 
example, using strong potentials one can get arbitrarily close to 100% efficiency. 



4 




Figure 2: Graphical representation of the reversible ratchet described in the text: the potential 
depends on two parameters, V\ and V2, which are the height of two barriers/wells (left) and 
they change along the path depicted on the center (V being the maximum height/depth of the 
barriers/wells). On the right, the shape of the potential at the four labelled points. 

To check the validity of the above theory and stress the differences between reversible and 
irreversible ratchets, we have studied in detail one example of each class. 

As an example of irreversible ratchet, consider the modulation of the potential in fig. |l| (left), 
i.e., V(x;t) = cos 2 {-Kt/T)V(x) with V(x) given by eq. (0). In this case, </>o is zero and the above 
theory cannot be applied. We have numerically integrated the Smoluchowski equation, eq. (||), 
using an implicit scheme with At = 10~ 5 , Ax = 0.002, 0.005, and the Richardson extrapolation 
method to correct inaccuracies coming from the finite Ax. The efficiency has been obtained using 
eqs. (§, (§), (0) and the results, as a function of F and for different values of T, are plotted in 
fig. H (left). The efficiency is maximum for T around 0.5 and it goes to zero as T increases. The 
maximum efficiency found by numerical integration is of the same order of magnitude as the one 
found for the randomly flashing ratchet. Notice, however, that we cannot explore with numerical 
experiments the whole space of parameters. 




Figure 3: Left (irreversible ratchet): numerical results for the efficiency of the ratchet consisting 
of the potential in fig. |l| (left) modulated by z(t) = cos 2 (7rf/T) as a function of the external force 
F and for different values of the period T: T= 0.00125 (O), 0.025 (□), 0.05 (O), 0.25 (x), and 0.5 
(A). 

Right (reversible ratchet): numerical and analytical results for the efficiency of the ratchet described 
111 fig. I for V = 5 a = 0.2 as a function of F and for different values of the period T: T= 1 (x), 
2 (O), 10 (□), 40 (O)- The thick solid line is the analytical result given by eq. ( pi] ) in the limit 
T — > 00 and F — > 0. Note that 77 is an increasing function of T in the reversible ratchet (right) as 
opposite to the irreversible case (left). 

On the other hand, let us consider the reversible ratchet represented in fig. ||. Here the 
potential depends on two parameters, Vi and V2, which are the heights/depths of two triangular 
barriers/wells of width a. The ratchet consists of modifying at constant velocity the parameters 
V\ and V2 along the path depicted in the same figure (center). This example is a modification of 



the one presented in ref . 1 1 . Now </> does not vanish and the above theory gives us the efficiency 



in the limit T — > 00 and F — > 0. For instance, for V = 5 and a = 0.2, we obtain O = 0.825, 
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/I = 0.094 and b — 3.74. The efficiency given by eq. ( |l~l|) is plotted in fig. || (right) and it is 
compared with the numerical integration of the Smoluchowski equation for different values of T. 
Notice the differences with the irreversible ratchet. Here the efficiency is an increasing function 
of T. The maximum efficiency, for the parameters corresponding to fig. |], is 26% and is almost 
reached for T — 40. The efficiency of this ratchet can be arbitrarily close to f 00% if one increases 
V. The reason is that the average mobility decreases exponentially with V, but the coefficient b 
and the integrated flow tfo remain finite. For instance, for V — 20 and a = 0.4, <fio = 0.999988, 
b = 6.89 and p, < 10 -7 , yielding a maximum efficiency of 99.85%. 

To summarize, we have calculated the efficiency of a randomly flashing ratchet with an asym- 
metric sawtooth potential. In order to find more efficient Brownian motors, we have also calculated 
the efficiency of deterministically driven ratchets, finding that the efficiency of reversible ratchets is 
much higher than the efficiency of irreversible ratchets. It is remarkable that the class of reversible 
ratchets involves potentials depending on two or more parameters O] and they differ from the 
models considered to date in the literature. Here we have shown that this new and non trivial 
class of ratchets is a real breakthrough regarding efficiency. 
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Tecnica (DGICYT) (Spain) Project No. PB94-0265. J.M.B. and F.J.C. acknowledge financial 
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